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I .  INTRODUCTION 


Spacecraft  which  contain  spinning  rotors  and  nominally  despun 
platforms  are  called  "dual-spin"  spacecraft.  In  such  a  dual-spin 
configuration,  the  rotor  provides  gyroscopic  stiffness  for  stability 
while  the  despun  platform  provides  an  oriented  platform  which  usually 
contains  scientific  instruments,  antennas,  solar  panels,  and  other 
components  which  must  be  oriented  in  a  "fixed"  direction. 

A  classical  rigid-body  analysis  of  the  attitude  stability  of  a 
spinning  single  body  spacecraft  indicates  that  a  state  of  steady  spin  of 
the  spacecraft  is  stable  if  rotation  is  about  its  principal  axis  of 
either  minor  or  major  moment  of  inertia.  This  classic  stability 
criterion  has  been  known  to  be  inadequate  since  1958  when  the 
unanticipated  instability  of  motion  about  the  minor  spin  axis  of  the 
first  U.S.  satellite,  Explorer  I,  was  observed.  The  explanation  of  this 
instability  is  credited  to  R.  N.  Bracewell  and  0.  K.  Garriott,  who 
modeled  that  spacecraft  as  a  semirigid  body  which  dissipated  energy,  and 
which  therefore  would  approach  the  state  of  minimum  energy  which 
corresponded  to  rotation  about  its  axis  of  maximum  moment  of  inertia.* 

In  1964,  Landon  and  Stewart  demonstrated  that  the  motion  of  the 
spin  axis  of  a  symmetric  dual-spin  spacecraft  with  a  despun  platform  may 
be  stable  if  the  spin  axis  is  the  axis  of  either  maximum  or  minimum 
moments  of  inertia  provided  an  energy  dissipation  device  (or  damper)  is 
placed  on  the  despun  platform.2  Independent  development  of  the 


2 


dual-spin  stabilization  concept  by  A.  J.  Iorillo  of  Hughes  Aircraft 
Company  extended  the  scope  of  analysis  to  axisymmetric  dual-spin 
vehicles  with  energy  dissipating  dampers  on  both  the  rotor  and 
platform.3  This  extension  allowed  for  the  analysis  of  spacecraft  of 
realistic  complexity,  since,  in  reality,  internal  motion  on  both  the 
rotor  and  platform  cause  energy  dissipation. 

Although  the  problem  of  attitude  stability  of  dual-spin  spacecraft 
has  been  studied  for  more  than  two  decades,  it  continues  to  be  an  area 
in  which  more  can  be  learned.  Several  different  methods  can  be  used  to 
study  the  attitude  motion  of  a  specific  spacecraft,  but  if  the  type  of 
configuration  permits,  a  conventional  linear  stability  analysis  is  a 
convenient  first  step.  A  linear  stability  analysis  of  a  specific  dual¬ 
spin  spacecraft  which  contains  arbitrarily  located  spherical  dampers  has 
been  presented  by  Laskin,  Sirlin,  and  Likins.4  They  also  obtained  cor¬ 
roborative  results  using  an  energy-sink  method.  In  applying  both 
methods,  small  nutation  angles  were  assumed.  An  alternative  approach  to 
the  energy-sink  method  is  available,  since  when  the  energy  dissipation 
is  slow  enough  to  justify  the  use  of  the  energy-sink  method,  a  pertur¬ 
bation  method  may  be  often  used  effectively.  An  analysis  of  this  type 
was  used  by  Cochran  and  Shu  to  analyze  the  nutational  motion  of  a  dif¬ 
ferent  spacecraft  configuration  and  they  obtained  very  good  results.5 

The  model  considered  in  Reference  4  is  of  the  "ideal"  type  for 
analysis  because  it  is  axisymmetric  and  the  energy  dissipating  devices 
on  the  rotor  and  platform  are  such  that  energy  can  be  dissipated 
without  changes  in  the  inertia  properties  of  the  spacecraft.  However, 
because  the  spacecraft  attitude  motion  is  nonlinear,  the  linear  analysis 
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does  not,  and  cannot  be  expected  to,  predict  correctly  the  attitude 
stability  of  the  spacecraft  with  respect  to  substantial  perturbations  in 
initial  conditions.  Laskin,  ^t_  ^1 .  ,  also  found,  via  simulation  on  a 
digital  computer,  that  in  some  cases,  stability  is  dependent  upon  the 
ratio  of  the  spacecraft's  axial  and  transverse  moments  of  inertia,  a 
parameter  which  does  not  appear  in  the  linear  analysis  equations.4 

This  thesis  presents  a  "perturbation  analysis"  of  the  attitude 
motion  of  the  type  of  spacecraft  investigated  in  Reference  4.  Following 
to  some  extent  Cochran  and  Shu,5  the  motions  of  the  spherical  bodies  are 
treated  as  perturbations  of  the  attitude  motion  and  the  generalized 
method  of  averaging  is  used  to  find  an  approximate  analytical  expression 
for  the  nutation  angle  for  the  case  of  a  constant-speed  rotor.  This 
expression  approximates  the  nonlinear  characteristics  of  the  spacecraft 
motion  and  also  contains  the  inertia  ratio  that  is  missing  from  the 
linear  stability  analysis.  The  expression  obtained  for  the  nutation 
angle  is  used  to  analyze  the  stability  of  the  spacecraft's  attitude 
motion. 

The  full  nonlinear  equations  of  motion  are  developed  in  Appendix  A 
and  numerical  solutions  of  these  equations  and  are  used  to  verify  the 
results  obtained  using  the  approximate  solution. 

The  equations  used  in  the  generalized  method  of  averaging  are 
developed  by  treating  the  changes  in  the  angular  moment  of  the 
spherical  dampers  as  perturbing  torques.  The  plat  form- fixed  components 
of  the  total  angular  momentum  of  the  spacecraft  are  used  as  dependent 
variables  along  with  the  components  of  angular  momenta  of  the 
dampers  and  the  rotor's  spin  axis  angular  momentum.  Euler  angles  which 


define  the  attitude  of  the  platform  are  then  used  to  make  a  change  of 
variables  which  transforms  the  angular  momentum  into  equations  of  the 
"normal  form"  required  to  use  the  generalized  method  of  averaging.6!7 
A  brief  description  of  the  generalized  method  of  averaging  is  given  in 
Appendix  B. 

The  generalized  method  of  averaging  is  used  to  produce  an 
expression  that  approximates  the  nutation  angle  of  the  spacecraft's 
attitude  motion  when  the  rotor  is  spinning  at  a  constant  speed  with 
respect  to  the  platform.  This  expression,  which  contains  all  the 
spacecraft's  inertia  and  damping  parameters,  is  used  to  determine 
regions  of  stable  and  unstable  spacecraft  attitude  motion  in  the  damping 
parameter  plane  for  several  spacecraft  inertia  parameters.  Results 
obtained  using  the  solution  for  the  "averaged"  nutation  angle  are  also 
compared  to  those  found  by  numerically  integrating  the  full  set  of 
nonlinear  equations  for  several  "constant-speed"  rotor  cases. 


II.  SPACECRAFT  MODEL 


The  type  of  spacecraft  considered  in  this  investigation  is  a  highly 
symmetric  dual-spin  configuration  such  as  that  depicted  in  Fig.  2.  The 
physical  model  of  this  type  spacecraft  consists  of  two  axisymmetric 
rigid  bodies,  a  nominally  despun  platform  and  a  rapidly  spinning  rotor, 
connected  by  a  shaft  running  along  their  collinear  axes  of  symmetry. 

Both  the  platform  and  rotor  contain  arbitrarily  located  spherical 
dampers.  These  dampers  may  be  thought  of  as  fully  loaded  fuel  tanks 
which  are  modeled  as  rigid  spheres  with  constant  surface  damping 
coefficients.4  Energy  dissipation  occurs  when  the  spherical  dampers 
rotate  with  respect  to  their  respective  containers  and  is  due  to  the 
torques  which  oppose  such  motion  of  the  dampers.  The  spherical  nature 
of  the  dampers  allows  for  energy  dissipation  with  no  change  in  the 
inertia  properties  of  the  spacecraft. 

In  Fig.  1,  the  Cx^x^x^  coordinate  system  has  its  origin  at  the 
center  of  mass  of  the  spacecraft  and  rotates  with  the  platform.  The 
unit  vector  triad  (e^e^.e^)  is  attached  to  Cx^x^x^.  The  angular 
velocity  of  the  rotor  with  respect  to  the  platform  is  £  *  fle^. 

In  this  analysis,  the  translational  motion  of  the  spacecraft  is  not 
considered  and  its  attitude  motion  is  assumed  torque  free. 
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III.  EQUATIONS  FOR  USE  IN  THE  GENERALIZED  METHOD  OF  AVERAGING 

Attitude  Motion  of  the  Spacecraft 
The  equations  for  the  attitude  motion  of  the  spacecraft  may  be 
developed  in  terms  of  _H,  the  angular  momentum  of  the  system.  The 
platform-fixed  components  of  JH  are 

H1  *  ItU)l  +  Ip(<*,pi+Wi)  +  IR^a)Ri+Tlii^  » 


«2  "  It“2  +  + 


H3  *  JPU3  +  JR^a>3+^  +  Ip(“p3+(,)3^  +  IR(f»p.j+ti>3+$)  * 

The  angular  momentum  of  the  dampers  has  platform-fixed  components 


Pi  *  Ip(“pi+Ui) 

,  for  i  *  1,2,3  , 

Ri  "  IR(uiRi+aii) 

,  for  i  31  1,2  , 

hR3  "  IR('J,R3+w3+^)* 

Hence,  Eqs.  (1)  may  be  rewritten  as 


H1  "  ItU)l  +  hPl  +  hRl> 


H2  *  It<02  +  hP2  +  hR2 


H3  *  JP“3  +  hP3  +  hR3  +  h  ’ 


[V*v> 


where  h  =  J 0(ui-+$),  is  Che  symmetry  axis  component  of  the  rotor's 
angular  momentum. 

If  the  no  angular  momentum  compor  'nts  of  the  platform  and  rotor 
dampers  are  combined  into  single  terms,  i.e., 

{hil  *  (hPi}  +  {hRil  ,  for  i  =  1,  2,  3  , 
then  the  following  expressions  for  ,  oi2 ,  and  result  from  Eq.  (3) 

uj1  -  (Hj-hj)/^  ,  ( 


u»2  -  (H2-h2)It 


u>3  »  (H3-h-h3)/Jp 


When  Euler's  moment  equation  is  applied  to  H,  the  following 
equation  is  obtained: 


M  *  H  +  ojxH 


Here,  H  is  the  time  rate  of  change  of  {h}  as  seen  in  the  Cx^x2x3 
system  and  is  the  external  torque  about  the  spacecraft's  center  of 
mass.  Since  the  motion  of  the  spacecraft  is  considered  torque  free, 
following  matrix  equation  may  be  obtained  from  Eq.  (6). 

{H}  -  -[2]  {H}  , 


where  {  h}  *  H2  &3]T  and 


0  -u>3  ui2 

[to]  ™  u>3  0 

-w2  w  ^  0 


-• ;  *  >  *>  *  *  v 


Equation  (7)  may  be  expanded  after  the  approximations  for  oi j ,  W2> 


and  103  are  substituted  to  obtain 


fil  - 

H2[“3^  -  T7>'h/Jp]  '  H2h3/JP  *  H3h2/It  ■ 

(9a) 

fi2  * 

'  T")  "  h/JP]  *  V3/Jp  -  Hjhi/i^ 

and 

Vint  -  W'c  ■ 

(9c) 

Letting 

A  -  I 

l  fi_  -  i_l  -  i_ 

3ljp  V  Jp  ' 

(10) 

Eqs.  (7)  may 

be  rewritten  as 

h  " 

H2A  -  H2h3/Jp  +  H3h2/Ic, 

(11a) 

*2  * 

-HjA  +  -  H3h1/Ic 

(lib) 

and 

h  * 

H2h1/It  -  Hjh2/Ie. 

(11c) 

Now,  refer  to  Fig.  2  which  shows  the  fixed  reference  frame  CXYZ 
with  the  constant  angular  momentum  vector,  JH,  aligned  with  the  Z  axis. 
The  spacecraft  axis  system  Cx^XjX^  is  obtained  from  the  CXYZ  system  by 
a  sequence  of  rotations  through  the  Euler  angles  Y ,  0,  and  t.  It  can  be 
seen  from  Fig.  2  that  the  Euler  angles  can  be  used  to  write  the 
components  of  H  in  the  form 


■  H  sin©  sin$, 
H2  *  H  sin©  cost 

and 

H„  -  H  cos© . 


(12a) 

(12b) 

(12c) 


By  differentiating  the  above  equations  with  respect  to  time,  one  finds 


■  H0  cos©  sint  +  sin©  cost, 


-  H0  cos0  cost  -  Hf  sin©  sint 


=  -H0  s  in0 . 

Equations  (11),  (12)  and  (13)  may  be  used  to  obtain  the  following 
equations : 


0  cos0  sint  +  t  sin©  cos® 


.  ^3 

A  sin©  cos<6  -  — —  sin©  cost  +  -=■—  cos0  , 
JP  Xt 


0  cos©  cost  -  t  sin©  sint  * 


^1  h2 

-A  sin©  sint  +  —  cos0  sint-  —  cos0 
Lt  Lt 


«  ^2  .  ^1 
0  sin©  *  —  sin©  sint  -  y—  sin©  cost. 

Lt  It 


Equations  (14b)  and  (14c)  may  then  be  solved  for  0  and  t,  respectively 


Hence, 


.  ^2  ^1 
0  *  —  sint  -  -j—  cost 


■A  -  t~  +  (h. sint+h„cost )/(!.  tan0). 


(18c 


13 


¥ 


H1  H2 

(—  sin4  +  ~  cos<t)/sin0 


By  substituting  values  for  ,  H2, 
attitude  equations 

0*0  +  ® (a) , 

4  ■  A  +  (5^(a) 

and 

¥  -  |-  +  &(a) . 


and  H3  from  Eqs .  (12),  one  gets  the 


(19a 

(19b 

(19c 


In  the  application  of  the  generalized  method  of  averaging  the  angles  $ 
and  ¥  are  referred  to  as  "fast”  variables  because  the  zeroth-order 
approximations  to  their  time  rates  of  change  are  non-zero.  The  angle  0 
is  called  a  "slow"  variable,  because  to  zeroth-order  it  is  constant.7 
If  the  terras  of  (a)  are  disregarded  and  X  »  —  ,  then  the  zeroth- 
order  approximation  for  the  Euler  angles  are 


0 


0o’ 

At  +  4. 


and 


(20a 

(20b 


¥  =»  Xt  +  ¥q.  (20c 

These  approximations  will  be  used  later  in  the  development  of  the 
zeroth-order  approximations  for  the  attitude  motion  of  the  dampers. 


Attitude  Motion  of  the  PlatformDamper 
To  develop  the  desired  equations  for  the  platform  damper  attitude 


motion  in  the  desired  form,  it  is  necessary  to  use  the  inertial 
components  of  angular  momentum.  If  { kp}  is  used  to  denote  the 
matrix  of  inertial  components  of  the  angular  momentum  vector  of  the 


V  >  J*' v*. . 


.  .  /.  _■-  _•.  \\\ «/, 
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platform  damper  and  a  3-1-3  rotation  through  the  angles  ¥ ,  0,  and  $ 
is  used  to  transform  from  inertial  components  to  platform-fixed 
components,  then 

{\}  "  [A]  {kp}  (21) 

where  [a]  is  the  transformation  matrix. 


cost  sint  0" 

*1  0  O' 

cos?  sin1?  0" 

[A]  3 

-sint  cost  0 

0  cos0  sin0 

-sin*?  cos?  0 

.  0  0  1. 

.0  -sinG  cos©. 

.  0  0  1. 

In  its  expanded  form  the  [ a] 

matrix  is 

(22) 


C08t 

-sint 

cos? 

COS0 

sin? 

cost 

+sint 

sin? 

cosG 

cos? 

sint 

sin© 

-sint 

cos? 

-sint 

sin? 

cost 

sin© 

-cost 

COS0 

sin? 

♦cost 

cos© 

cos? 

sin0 

sin? 

-sin© 

cos? 

cosG 

(23) 


Euler's  moment  equation  may  be  applied  using  the  angular  momentum  of  the 
platform  damper  to  obtain  the  vector  equation, 


h  +  wx  h  3  T  =  -Coj 
-P  -  ~P  -P  P“P 

By  writing  Eq.  (24)  in  matrix  form  and  substituting  { ojp } 
one  may  obtain  the  results,  ✓ 


(24) 

f-  (V  -M  ■ 

P 


(25) 


If  Eq.  (21)  is  now  differentiated  with  respect  to  time,  a  second 
equation  for  fh^}  is  found  to  be 


ty  -  [a]  {y  *  [a]  {y. 

However,  it  can  be  shown8  that 
[A]  »  >  M  [A]. 

Hence,  Eq.  (26)  may  be  written  in  the  form, 

ty  --[«]  [A]iy  +  [ a]  {y. 

Mow,  by  equating  Eq.  (28)  and  Eq.  (25)  and  using  Eq.  (21)  for 
{hp},  one  may  obtain  an  expression  for  {kp};  i.e., 

IM  • 

A  solution  to  Eq.  (29)  may  be  assumed  in  the  form, 


-0pt/Ip 


(M  ■  (“pi  e 


If  Eq.  (30)  is  differentiated  with  respect  to  time,  a  second  form 
for  {kp} ,  viz. , 

ty  -  {ye'CpC/Ip  -i{kpj. 

Mow,  by  equating  Eq.  (29)  and  Eq.  (31),  one  finds  that 

.  Cpt/l_  T 

(ipl  -  V  P  P[A]T(»!  . 

Equation  (32)  may  be  expanded  to  obtain 
Cpt/Ip 

opi  -  Cpe  (auu,1  +  a21u,2  +  a31  0,3), 

CPC/IP 

°P2  "  CP*  (a12W  1  +  a22u,2  +  a32  a,3) 
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Cpc/Ip 

Opa  =  Cpe  ( ai 3^ i  +  a23a,2  +  a33  “3)  (33c) 

where  a^  j  are  Che  elements  of  the  [A]  matrix. 

If  the  approximations  for  0,  $,  and  given  in  Eqs.  (20)  and  the 
approximations  for  uii,  ui2,  and  013  given  by  Eqs.  (17)  are  used  in  Eqs. 
(33),  it  can  be  shown,  after  some  algebra,  that 

dpi  *  Cp  A  z  e  sinf  ,  (34a) 


C  t/I 

oip2  *  ~cp  A  z  e  cosf  (34b) 

and 

CPC/IP 

(Xpa  *  Cp[Xz2  +  x(H3-h)/Jple  ,  (34c) 


where  z  *  sin0  and  x  *  cosQ  . 


The  solutions  to  the  above  equations,  integrated  from  to*0  to  ti*t. 


are 

CpAz  C-t/Ip  Cp 

Op!  *  -5 -  (e  (f“  sin*  -Acos'f)  +  X}  ♦  apl0,  (35a) 

(r)2**2  P 

P 


-C-Az  Cpt/Ip  cp 

ip 2  *  — -  {e  ( cosf  +X  s  inf  )  -  — }  +  a 


1P 


P20 


(35b) 


and 


P3  P 


Cpt/I 

Ip{Xz2  ♦  xH,-h)/Jp}(e  -1)  +  ap30. 


(35c) 


When  used  in  Eq.  (30),  the  above  expressions  give  the  inertial 
components  of  the  platform  angular  momentum  as  functions  of  time  and  the 
fast  variable  '{'■Xt.  These  components  are 
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cpAz  c  -C  t/r  -c  t/i 

kpi  -  -p- -  (j-  sin*  -Xcos*  +Xe  )  +  apl()e  , 

P>2.,o  P  (36a) 


(rr-*2 


and 


^P2 


kp3 


>CpAz  Cp 


r  f—  cos*  +Xsin*  -  —  e 

It 1)2+x2  p  p 


Cp  -Cpt/Ip  -Cpt/Ip 

)  +  aP20e 


-C  t/1  -cpt/i 

Ip[Xz2  +  X(H3-h)/Jp](l-e  P  P)  +  ap30e  P  P. 


(36b) 


(36c) 


Armed  with  an  expression  for  Che  inerdal  components  of  angular 
momentum,  one  may  find  the  relative  components  from  Eq.  (21),  which  in 
component  form  are 

hPl  *  allkPl  +  al2kP2  +  a13kP3*  (37A) 

hP2  *  ^l^Pl  +  a22kP2  *  ^^PS  (37b) 

and 

hP3  *  *31^1  +  a32kp2  +  *33^3  *  (37c) 

If  Eqs.  (37)  are  expanded  using  the  approximations  for  0 ,  * ,  and  $ 
in  the  transformation  matrix  and  the  transient  terms  in  {kp}  are 
neglected,  then  it  can  be  shown,  again  with  some  algebra,  that: 

CpZA  Cp 

hpi  ‘  Tcprpr»?  t-Xc°st  -  t;  * sintl 


+  Ipz  sin*  [Xz2+x(H3~h)/Jp]  , 


CpZA  Cp 

hp2  ‘  (cp/i7^  •  t;  * c°,tl 

+  Ip  zco8$  (Xz2  +  x(H3~h)/Jp] 


(38a) 


(38b) 


and 
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^2  ^2. 

hp3  =  (cpyip)2"^"7  *  1 P  x[Xz2  +  *(H3-h)/Jpl.  (38c) 

The  above  expressions  give  Che  angular  momentum  components  of  the 
platform  dampers  as  functions  of  the  slow  variable  0  found  in  x  and  z 
and  the  fast  variables  *  and  ¥. 


Motion  of  the  Rotor  Pampers 

The  equations  for  the  rotor  damper  motion  may  be  developed  in  the 
same  manner  as  those  for  the  platform.  If  {kR}  denotes  the  inertial 
angular  momentum,  then  the  same  transformation  exists  for  the  rotor 
dampers,  i.e., 

< 

{hR}  *  [A]  {kj  •  (39) 

If  a  vector  {$}  is  defined  such  that 

{$}  -  (0,  0,  $)T  ,  (40) 


then  Eqs.  (2b)  and  (2c)  can  be  written  in  matrix  form  as 

tv  "  M“iJ  +  xrW  +  xrW* 


(41) 


Euler's  moment  equation  applied  to  the  angular  momentum  of  the  rotor 
gives 


-R  +  ^R 


(42) 


Switching  to  the  matrix  notation  and  substituting 

{ti)R}  ■  y-  {hR}  -  fu}  -  {$},  into  Eq.  (42)  yield 
R 

{\}  -  -VXR  { hR 1  +  CRW  +  Cr{  '  I«1  {\}  •  <«) 


By  differentiating  Eq.  (39)  with  respect  to  time>one  obtains 


{\\  =  [A]  {kj  *  [A]  {kR}  =  -  [U]  [A]  {kR}  +  [A]  {kR}  ( 

Now,  by  equating  Eqs.  (43)  and  (44)* it  follows  that 

{\}  ■  '  V!R  WJ  +  CR[A]T  M  +  CR[A]T  {♦}■  ( 

As  with  the  platform  dampers,  the  solution  for  {kg}  is  assumed  to  be 


-V^R 


("rI  ■  W8 


By  differentiating  Eq.  (46)  with  respect  to  time  and  equating  the 

result  to  Eq.  (45), one  may  obtain  the  following  expression  for  {a  }; 

K 

CUC/lO  T  C_t/I_  _ 

{c^}  »  CRe  R  R[A]T  {oa}  +  CRe  R  R[A]T  {$}  .  (47) 

If  the  approximate  expression  for  |u>}  and  the  approximations  to  0 ,  , 

and  ¥  are  used  in  the  transformation  matrix  in  Eq.  (47)  then  the  results 


for  the  a_ .  are 


•  V/XR 

aRj  *  CRz(A+$)e  sin¥. 


;  .  V71* 

R3  *  -CRz(A+$)e  cos¥ 


aR3  *  CR [X z^  +  x(H,-h)/JB  +  x  |]e 


CRC/IR 


3 


(48a) 


(48b) 


(48c) 


Equations  (48)  may  be  integrated  from  tQ*0  to  tj*t  to  get  the 
following  results: 

C  «(A++)  C  t/T  cn 

Ri "  rcjry-  u  te  9in'i' _x  c°8f)  +  +  “rio* 


(49a) 


-Crz(A.»)  CRt/I  CR 


aR2  "  (C_/I  K  ♦  t 8 


\  Xs 


cos¥  +Xs 


'  rt  88r2o  ..... 
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and 


.  CRC/lR 

aR3  *  rR  ^ Xz 2  +  x  (H3_h)Jp  +  x  <H(e  -l)  +  aR3Q  (49c) 

When  used  in  Eq.  (46),  Eqs.  (49)  provide  the  inertial  components  of  the 
angular  momentum  of  the  rotor  damper.  These  are 


C.z(A+f)  C 

kRi  •  ic~/iRy  t  1-^ 


]  ♦  « 


-vn* 


RIO 


(50a) 


and 


CRz(A  <j>)  CR  .  CR  Cr^r  CR20 

^2  *  rc'Trwv  cos*  Xslnf  -  f- e  +  «R2oe 

™  CW  A  rR  \  RZ0  (50b) 


.  -CRt/l  -C  t/1 

kR3  *  LR  ^Xz2  +  x^H-i-h^Jp  +  ](  1-e  )  +  ar„„e 


R30 


(50c) 


From  these  expressions  for  the  inertial  angular  momentum  components,  one 
may  find  the  platform-fixed  components  of  the  angular  momentum  by  using 
Eq.  (39)  with  the  zeroth-order  approximation  to  the  transformation 
matrix  and  {<u}  .  The  results  are 


— C^z(A*h^)  CR 

hRi  •  (c JiJi  *  x*  tr x  !in*  * XcostI 

K  K  K 

♦  Ip  z  sinAtfAz2  +  x  (H^-hJ/Jp  +  x  $  ]  ,  (51a) 


CRz(A+$)  Cr 

hR2  •  (c./i.W  *  57  (x,in* '  r  *  cos*1 

R  R  R 

+  Ip  z  cosAt[X^2+x(Hp-h)/Jp  ♦  x  $]  (51b) 

and 

c2*2(A+i)/l 

hR3  *  TcTTTT?  *  &  +  rRx[Xz2  +  x( ( H3-h ) / Jp  ♦  x  ♦]  (51c) 

R  R 


.v  > 


i  Mif  fr~i  ^  1 


Equations  (51)  are  zeroth-order  approximations  to  components  of  the 


angular  momentum  of  the  rotor  damper  as  functions  of  the  slow  variable  0 


and  the  fast  variables  $  and  'll  . 


Transformation  of  the  Equations  to  Normal  Form 


To  apply  the  generalized  method  of  averaging,  it  is  necessary  to 


have  equations  of  the  form 


{*}  3  e{Xj}  ♦  £2{X2}  +  ... 


(52a) 


{v}  *  {Yq}  ♦  +  e2{Y2}  + 


(52b) 


where  e  is  a  small  constant,  { x}  is  a  vector  of  slow  variables,  { y}  is  a 
vector  of  fast  variables,  and  the  vector  functions  { X £ }  and  {Yj}  are 
periodic  functions  in  the  elements  of  { y}  with  a  period  of  2ir  ? 


To  arrive  at  the  desired  normal  form,  one  must  combine  the 


angular  momentum  solutions  to  get 


ij  =■  -  A  cost  [- 


C  z(A+$) 


VV  + 


!„)*  +  AZJ 


c|zA/Ip  C|z(a4)/Ir  , 
-  x  sin$  [(Cp/Ip)z  \\l  +  (Cr/Ir)2  T  \T1 


+  z  sin$  [A z2  +  x(H2-h)/JR](lp+IR)  +  IRxz$  sin$ ,  (53a) 


CpzA  Crz(A+$) 

h-2  ■  Xsin<  [(c^ipy  .  j?  + 


.  ,  .  cI'(a**)/ir, 

co8*  ITcJTpz  7  x'i  *  Ccji-T-nji 


+  z  cos$  (Az2  +  x(H_-h)/J  ](l  +I_)  +  I  xz$  cos<t  (53b) 

j  r  r  K  K 


CpZA/Ip 

rc;^— 


c|z(a^)/ir 

“V1^  *  : 


+_X?J 


+  xtXz2  +  X(H3-h)/jp](lp+  IR)  +  IR  x2* 


(53c) 


The  above  equations  for  h^  and  h2  along  with  the  approximation  for  $ 
given  by  Eq.  (15a)  may  be  used  to  obtain  the  result, 


g  _  X  r 

It  KCp/Ip)2  + 


CRz(A+$  ) 

w  * xJ 


Now,  it  is  assumed  that  the  inertias  of  the  dampers  are  equal  and 
a  ■  lp/It  *  IR/It  adopted  as  the  small  parameter.  Since 


CP/Tt  *  (lP/It)(CP/Ip)  ~a  CP;iP 


CR/It  *  (lR/It)(CR/lR)  =*  a  CR/IR, 


Eq.  (54)  may  be  written  in  the  form 


CpXzA/Ip 

0  •«  (ryi^-T 


CpXzfAH'/Ip 

F"  *  Tyyr  Fid  • 


Now,  because  z  and  A  in  the  above  equations  are  functions  of  x,  it  is 
chosen  as  the  slow  variable.  Prom  x  *  cos©  it  follows  that 


x  *  -0  sin©  =  -©  z 


I 


Hence,  0  =»  -  x/ z  and  Eq.  (57)  may  be  replaced  by 


Equation  (59)  is  one  of  the  desired  equations  in  normal  form.  Th 
other  two  are 

I  =  A  +  O'  (a)  (60 

and 

V  =  X  +  a)  (60 

Equations  (59)  and  (60)  are  in  the  required  normal  form  for  the 
generalized  method  of  averaging.  Only  Eq.  (59)  is  expanded  because  it 
is  the  only  one  required  to  investigate  nutational  motion. 


IV.  FIRST  ORDER  SOLUTION  FOR  THE  NUTATION  ANGLE 


The  first-order  solution  for  the  nutation  angle,  in  terms  of  x,  c 
now  be  developed  by  applying  the  generalized  method  of  averaging  to  Eq 
(59).  Initially,  the  following  nondimens ional  ratios  must  be  defined: 


WP  '  <5£  +  X2 


MR  62  +  A2  ’ 

Then  Eq.  (59)  may  be  written  as: 

x  »  ot{-ypz2A  -  yRz2(A+|)] 

At  this  point  one  may  apply  the  generalized  method  of  averaging  to 
Eq.  (62).  Referring  to  the  procedure  in  Appendix  B, 

,  2tt  2ir 


1  All  All  _  _ 

Ax(x)  =  /  /  X^x,*  ,¥)d*d¥ 


0  0 


where 


X^(x,t,V)  *  -TJpz2A  -  p  pZ2  (A+$ ) 


Since  X^  does  not  contain  periodic  terms  in  <t  and  ¥ ,  the  averaged 
equation  to  first  order  in  the  new  parameter  variable  o  is 

x  *  a[~iipZ2  A  -  u  ?(A+$)]  (( 

where  z  and  A  are  obtained  by  replacing  x  with  x  in  z  and  A. 

Equation  (65)  may  be  solved  exactly  by  separating  variables  and  using 


( 


b  -  —  • 

Now,  by  separating  variables,  using  partial  fractions  and  integrating 
Eq.  (71)  with  x  *  xq  at  t  ■  0  the  following  solution  may  be  obtained. 

t  *  -C.  £n[  (l-x)/(l-x  )]  +  C«  Jln[(l+x)7(l+x  )] 

1  o  z  o 

+  C,  £n[(a+bx)/(a+bx  )]  ,  ( 

j  o 

where 

C1  "  2(a+b)  * 

r  -  1 

C2  2(a-b) 


V.  STABILITY  CONDITIONS 


If  the  spacecraft's  nominal  attitude  motion  i3  asymptotically 
stable,  the  "averaged"  nutation  angle  will  decay  to  zero  as  time 
increases.  Considering  Eq.  (71)  for  asymptotic  stability,  one  sees 
that  it  is  necessary  that  x  should  be  positive.  This  criteria  may  be 
expressed  as 

(l-x2)(a+bx)  >  0.  (75 


Since  the  term  (1-x2)  is  always  positive  for  nutation  angles  greater 
than  zero,  it  is  necessary  that 


a  +  bx  >  0  . 


The  substitution  of  the  explicit  expressions  for  a  and  b  into  Eq.  (76) 
provides  the  inequality, 

|(uRp-upa) 

x  >  H/It(up+uR)(o+p-l)  *  (7? 

If  the  averaged  value  of  x  is  used,  then  *  Hx,  and  the  condition  for 
stability  is 

H3/lc(upmR)(o+p-l)  >  $(uRp-Vpo)  •  (78 


If  the  stability  of  the  nominal  spacecraft  motion  in  which  the  platform 
is  inertially  fixed  and  the  rotor  spins  at  a  constant  rate  is 
considered,  then,  for  small  nutation  angles 


H3/It  =  H/It  -  o*fl  , 


where  a*  *  a  +  a  and  J  *  !1  is  a  constant,  then  Eq.  (78)  may  be  put 
into  the  form, 

Up(o+p-l  +  o/o*)  +  u^Ca+p-l-p /a*)  >  0  . 

For  a  specific  set  of  inertia  parameters,  Eq.  (80)  can  be  used  to 
define  regions  of  stability  in  the  6p  -  6^  plane. 


(80) 


VI.  RESULTS 


Small  Nutation  Angles 

If  condition  (80)  is  examined  in  the  limiting  case  where  6g+0,  the 
requirement  for  stability  is  found  to  be 

a  +  p  +  a/a*  >  1  ,  (81) 

which  is  always  true  in  the  small  "a"  approximation.  This  case 
corresponds  to  energy  dissipation  on  the  platform  only  and  the  motion 
should  be  stable.2 

In  the  other  limiting  case  where  5p+-0,  the  requirement  for 
stability  is 

a  +  p  >  1  +  p/a*.  (82) 

A  careful  examination  of  condition  (82)  reveals  that  a  value  of  o>l 
would  insure  stability  of  nutational  motion  in  this  case.  This  condition 
corresponds  to  an  oblate  spacecraft  (i.e.,  one  that  spins  about  its  axis 
of  maximum  moment  of  inertia)  which,  according  to  previous  results,  also 
exhibits  stable  behavior.1 

When  energy  dissipation  is  present  on  both  the  rotor  and  platform, 
condition  (80)  can  be  examined  for  specific  inertia  parameters.  Figures 
3  through  7  are  stability  diagrams  for  6 p/Q  vs.  <5g/£2  for  different 
spacecraft  parameters.  Fig.  3  shows  a  critical  case  where  o  and  p  both 
equal  0.5. 

Some  qualitative  results  can  be  obtained  by  examination  of  Figs.  4 
through  7.  These  figures  indicate  that  if  the  rotor  inertia  parameter 
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Figure  3.  Stability  plot  in  the  damping  parameter  plane  for 
ct-0.05,  p-0.5,  a-0.5. 
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is  larger  than  the  platform  inertia  parameter,  then  more  of  the  area  in 
the  region  considered  predicts  stable  motion,  while  if  the  platform 


inertia  parameter  is  larger,  then  more  of  the  region  predicts  unstable 
motion.  In  Figs.  4  and  5,  the  platform  ratio  is  the  larger  and  most  of 
the  area  predicts  unstable  motion,  especially  in  Fig.  5  where  the  rotor 
parameter  is  less  than  the  critical  value  of  0.5.  In  Figs.  6  and  7,  the 
rotor  ratio  is  the  larger  of  the  two  and  most  of  the  area  predicts 
stable  motion,  especially  in  Fig.  7  where  the  rotor  parameter  is  larger 
than  0.5.  It  should  be  noted  here  that  Figs.  3  through  7  are  very 
similar  to  the  stability  diagrams  found  by  using  a  linear  analysis.4 
However,  the  results  of  this  analysis  contain  the  effects  of  inertia 
parameters  for  both  the  rotor  and  platform,  while  the  linear  analysis 
results  depend  only  on  the  rotor  inertia  parameter. 

Nonlinear  Stability 


Another  result  of  the  generalized  method  of  averaging  is  that  an 
analytical  approximation  to  the  nutational  motion  can  be  found  from 
Eq.  (74).  Figures  8  through  11  show  how  the  approximate  solution  of  Eq. 
(74)  compares  to  the  "exact"  numerical  solution  of  the  full  nonlinear 
equations  (see  Appendix  A).  The  full  nonlinear  equations  were  solved 
using  a  fourth-order  Runge-Kutta  method. 

To  compare  the  results  of  the  approximate  solution  to  the  Runge- 
Kutta  solution,  two  cases  are  considered;  one  with  the  platform  inertia 
parameter  larger  than  the  rotor  inertia  parameter  and  one  with  the  rotor 
inertia  parameter  the  larger.  The  predicted  subcases,  one  stable  and 
one  unstable,  of  each  case  are  considered.  The  nonlinear  equations  are 


NUTATION  ANGLE  (DEG) 

0.00  10.00  20.00  30.00  40.00 
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integrated  with  $  a  0  and  initial  conditions  of  =  0.2  rad/sec, 

(2  *  1.0  rad/sec  and  all  others  zero. 

Figures  8  and  9  show  the  results  obtained  when  the  platform  inertia 
parameter  is  larger  than  the  rotor  inertia  parameter  and  the  approximate 
results  are  in  good  agreement  with  the  numerical  results.  Figures  10 
and  11  show  the  results  for  the  case  when  the  rotor  parameter  is  larger 
than  the  platform  inertia  ratio  and  again  the  results  agree  well, 
especially  in  the  stable  case  of  Fig.  10.  In  Fig.  11,  the  transient 
nonlinearities  seem  to  increase  the  difference  between  the  approximate 
and  numerical  results.  In  all  cases,  the  difference  between  the  two 
solutions  appears  to  be  a  linear  function  of  time. 

Figures  12  and  13  show  the  results  obtained  for  cases  of  no  damping 
on  the  rotor  (Fig.  12)  and  no  damping  on  the  platform  (Fig.  13), 
respectively.  These  figures  show  the  effects  of  neglecting  the  transient 
terms  in  Eqs.  (36)  and  (50).  When  there  is  no  damping,  the  oscillatory 
terras  never  disappear.  However,  that  these  terms  do  eventually  vanish  if 
only  a  small  amount  of  damping  is  present  is  illustrated  by  Fig.  14. 
the  transient  terms  die  out  very  rapidly  when  significantly  large 
damping  constants  are  used 

Figures  15  and  16  demonstrate  the  results  of  "fixing”  the  damper  on 
the  rotor  and  platform,  respectively.  These  simulation  results  were 
obtained  by  setting  the  damping  ratio  to  zero,  while  at  the  same  time 
setting  u>R1  -  -  ui^  -  0,  for  Fig.  15,  and  u>pi  -  u>p2  -  up3  *  0,  for 

Fig.  16,  in  the  full  set  of  equations.  Figures  15  and  16  are  analogous 
to  Figs.  12  and  13,  respectively.  Again,  the  analytical  and  numerical 
results  compare  quite  well,  especially  for  the  stable  case. 
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a-0.6,  6  =0.5,  5  *0.0. 

r  R 


*  <* 


NUTATION  ANGLE  (DEG) 


VII.  CONCLUSIONS 


The  attitude  motion  of  a  model  of  a  symmetric  dual-spin  spacecraft 
containing  spherical  dampers  on  both  its  rotor  and  platform  has  been 
investigated.  A  perturbation  method  which  treats  the  effects  of  the 
motion  of  the  dampers  as  perturbing  torques  has  been  used  in  conjunction 
with  the  generalized  method  of  averaging.  A  stability  analysis  was 
conducted  for  the  constant  speed  rotor  case  and  produced  results  that 
agree  well  with  the  linear  stability  analysis  described  in  Ref.  4. 
However,  the  analysis  of  this  thesis  is  more  general  in  that  more 
information  regarding  the  inertia  parameters  is  obtained  than  is 
possible  using  a  conventional  linear  analysis.  The  generalized  method 
of  averaging  also  provides  an  approximate  solution  for  the  nutational 
motion  which  is  not  restricted  to  small  angles.  This  solution  agrees 
well  with  results  obtained  by  numerically  integrating  the  full  set  of 
nonlinear  equations. 
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APPENDIX  A 


DERIVATION  OF  THE  SPACECRAFT  EQUATIONS  OF  MOTION 


APPENDIX  A 


The  equations  of  motion  of  the  spacecraft  are  derived  by  applying 
Euler's  moment  equation,  to  the  spacecraft  as  a  whole,  then  to  the  rotor 
and  platform  dampers  individually,  and  then  to  the  rotor  subsystem. 

The  angular  momentum  of  the  entire  system  about  its  center  of  mass 


+  Ipiaipj^+w^)  +  In^a,jQ+a1^  >  (A-la) 


H2  *  Itu2  +  Ip(a>p2+a>2)  +  Ir^i,R2+w2^  (A-lb) 


H3  *  JpU>3  +  JR(a>3+^)  ♦  Ip(Up3+m3)  +  IR(toR3+t»)3+|)  (A-lc) 


ut  i;  ■  h  *  h *  h  and  ■  jP  *  jR  *  h  *  h 

Then,  in  terms  of  Is*  and  Ic*,  one  has 


^“l  +  V»Pl  +  V*R1 


(A-2a) 


«2  “  Z*w2  +  V*P2  +  IRUJR2 


(A-2b) 


H3  “  ^“3  +  Vp3  +  Vs 3  +  (W* 


(A-2c) 


Now,  by  applying  Euler's  moment  equation  to  Eqs.  (A-2),  one  finds 


M  ■  H  +  o)xH 


(A-3) 
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Since  the  spacecraft  is  not  affected  by  external  torques 


o 

H  +  ojxH  =  0  .  (A-4) 

The  component  form  of  Eq.  (A-4)  is 

X*“l  +  Vpi  +  VrI  +  u2H3  ■  “3H2  =  0  (A"5a) 

Xt“2  +  IP“P2  +  IR.“r2  +  U)3H2  "  “iH3  =  0  (A-5b) 

Xs“3  +  XP^P3  +  rR“R3  +  (JR+V*  +  “lH2  '  «2H1  =  °  •  (A*5c) 


By  substituting  for  ,  and  and  simplifying  Eq.  (A-5)  one 

gets 

^“l  +  VP1  +  VR1  +  (JR+IR^“2  +  ^Xs”XtPM2W3 

+  Ip((Dp3w2  -  “p2“3^  +  IR^a)R3w2_a>82a>3)  “  °»  ^A_6a^ 

I*i2  *  Ipip2  +  Vr2  -  <VVK  +  (l*-I*)o)1a.3 

+  XP^,Pl<‘,3-uP3a,l)  +  ir(«ri“3“  "  0  (A-6b) 

and 

X^3  +  ^P3  +  IR“R3  +  (JR+IR)  ♦ 

+  Ip(up20)1-u)p1a)2)  +  IR(“R2u)i~  U)Rla,2^  *  °*  ^A_6c^ 

The  plat  form- fixed  components  of  angular  momentum  for  the  platform 
dampers  can  be  written  as 


1  m 

pi 

(A-7a) 

'P2 

IP(“p2*”2) 

(A-7b) 

‘P3  ” 

Ip(up3+  “3>. 

(A-7c) 

Application  of  Euler  s  moment  equation  to  the  platform  damper 


yields 


M  =  h  +  (oxhp  =  Tp 


(A-8) 


where  T_pis  a  3x1  matrix  of  the  components  of  the  torque  applied  to 
the  platform  dampers.  Equations  (A-8)  may  be  written  in  component  form, 


V^lV  +“2hP3  Ui3hP2  =  TP1 


Ip(ip2+m2)  ♦  o)3hpl  -  tl,1hp3  =  Tp 


XP (aJp3+<^3)  +  “lhP2  ~  “2hPl  =  TP3 


(A-9a) 


(A-9c) 


(A-9c) 


One  may  replace  Tp^  *  ^or  *  ~  1,2,3  and  substitute  for  hp^ 


for  i  ■  1,2,3  to  get  the  equations. 


IpCip^)  +  IP(a)P3a)2“wP2u,3)  +  CP“P1  *  °> 
Ip(ip2^2)  +  Ip(a>pla)3-<Up3a)1)  +  Cpa)p2  -  0 


XP(<^P3+<^3)  +  IP(a,P2“l'wPlU)2)  +  CP“P3  =  ° 


(A-10a) 


(A-lOb) 


(A-lOc) 


Similarly,  the  angular  momentum  for  the  rotor  dampers  can  be  written  as 


hRl  =  VU)R1+U,1)’ 


hR2  “  IR(o,R2+U)2) 


1R3  *  V“R3+ 


(A-lla) 


(A-llb) 


(A-llc) 


and  the  application  of  Euler's  moment  equation  to  the  rotor  dampers 
provides 


(A-12a) 
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V'W  *-“2hR3  •"3ha2  *  TR1  • 

IR^“R2+U2^  +  “3hRl  "  “l1^  “  TR2  (A-12b) 


and 


V“R3+  V 

♦>  '“lhR2 

ot 

x: 

IM 

3 

1 

3  T 

R3 

(A- 

■12c) 

By  using  TRi  » 

CR“Ri'  1  ‘ 

12  3 

and  substituting 

hRi  for 

i  =* 

1,2,3,  one  may 

put  Eqs  (A-12)  in  the  form. 

+  IR(l0R3UJ2"u 

R2U3)  + 

xrK  + 

cr“ri  * 

0, 

(A- 

■13a) 

f 

ro 

3 

2 

3 

& 

M 

+ 

RS^l^ 

V“l + 

CR“R2 

0 

(A- 

■13b) 

and 

V^RS+V^  +  V«R2“l 

"WR1U)2) 

+  CR“R3 

*  0. 

(A- 

•13c) 

By  considering  the  rotor  as  a  system,  one  may  derive  three 
equations  for  the  rotation  motion  of  the  rotor;  however,  only  the 

A 

equation  corresponding  to  the  rotation  about  e^  is  needed.  This 
equation  corresponds  to  the  tenth  degree  of  freedom  of  the  system. 

If  {R}  is  used  to  denote  the  3x1  matrix  of  platform-fixed 
components  of  the  angular  momentum  of  the  rotor  and  Jt  as  its  transverse 
moment  of  inertia,  the  following  equations  may  be  written: 


+  IP(WP14W1)  +  IR(uRl+tI,l>  * 

(A-14a) 

+  IP(“P2^2)  +  IR(uR2+<‘,2) 

(A-14b) 

and 


R3  *  JR(u3+$)  +  Xp (uj p3  +u) 3 >  +  . 


(A-14c) 


Euler  s  moment  equation  may  be  applied  to  Eqs.  (A-14)  to  get 


M  *  R  +  ojxR 


Considering  only  the  third  component  of  Eq.  (A-15),  one  gets 
(JR+lR)(“3+^  +  ¥r 3  +  IR^aJR2U>l~u>Rlai2)  “  0  ♦ 
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The  basic  idea  behind  the  method  of  averaging  is  to  derive  an 
approximate  solution  of  the  nonautonomous  system  by  considering  in  its 
place  an  associated  average  system,  which  is  autonomous.6  For  the 
generalized  method  of  averaging  to  be  used,  the  differential  equations 
to  be  solved  must  first  be  cast  in  "normal  form."6  The  algorithm 
described  below  applies  to  the  system  of  differential  equations  in 
this  thesis  which  may  be  expressed  as 


x  *  aXj (x,$ ,¥ ) , 


(B-la) 


*  ■  A(x)  +  aYj  (x,4  ,¥ ) 


CB-lb) 


and 


¥  -  X(x)  +  aZ^x,*  ,Y)  , 


(B-lc) 


where  Xj ,  ,  and  Z\  are  periodic  with  a  period  of  2ir  in  the  angles 

and  Y,  and  o  is  a  small  constant.  An  approximate  solution  to  the 
system  which  is  valid  through  first  order  in  a  may  be  obtained  by 
solving  the  "averaged"  equations: 


x  ■  a  A^ (x) , 


(B-2a) 


♦  ■  A(x)  +  a  Bj (x) 


(B-2b) 


and 


A(x)  +  a  (x) . 


(B-2c) 
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m 
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2tt  2tt 

/  /  X  (x,$  ,¥)d*df  , 

0  0 

2it  2tt 

/  /  Y  (x,$  ,Tf)d*d?  , 

0  0  1 

2ir  2ir 

J  /  Z.Cx.t.lOdtdY  , 

0  0  1 


A(x)  «*  A(x) 


A(x)  ■  A(x) . 

The  transformation  from  the  averaged  system  is 
x  a  x  +  a  Uj  (x,$  ,Y) , 


where 


$  *  *  +  a  (x,$ ,?) 


T  ■  ?  +  a  Wj  (x,$  ), 


.  2ir  2tt 

m^Z  !  I  (X  -A  )d*<ff  , 
A  w  o  0 

.  2it  2ir 

».  -  A  f  f  <Y.-B.)d*df 

1  **0  0  1  1 


.  2it  2ir 

w!  “  A  /  /  (Z  -C.)d*df  . 

*  w  0  0 


(B-3a 


(B-3b 


(B-3c! 


(B-4a' 


(B-4b: 


(B-5a' 


( B-5b' 


(B-5c: 


(B-6a' 


(B-6b 


(B-6c 


In  the  present  application,  the  expression  for  x  is  the  one  of  dominant 


interest . 


